On the relation between decoherence and spontaneous symmetry breaking 
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We have recently shown that there is a limit to quantum coherence in many-particle spin qubits due 
to spontaneous symmetry breaking. These results were derived for the Lieb-Mattis spin model. Here 
we will show that the underlying mechanism of decoherence in systems with spontaneous symmetry 
breaking is in fact more general. We present here a generic route to finding the decoherence time 
associated with spontaneous symmetry breaking in many particle qubits, and subsequently we apply 
this approach to two model systems, indicating how the continuous symmetries in these models are 
spontaneously broken and discussing the relation of this symmetry breaking to the thin spectrum. 
We then present in detail the calculations that lead to the limit to quantum coherence, which is due 
to energy shifts in the thin spectrum. 



INTRODUCTION 

Recently we have shown that spontaneous symmetry 
breaking imposes a fundamental Hmit for the time that a 
large spin system can stay quantum coherent. This uni- 
versal time-scale is ispon — 27TNh/{kBT), given in terms 
of the number of microscopic degrees of freedom TV, tem- 
perature T, and the constants of Planck (ft) and Boltz- 
mann (fcs) 1]. We analyzed this quantum decoherence 
process in terms of the exactly solvable Lieb-Mattis spin 
model, which is known to describe a symmetry broken 
macroscopic antiferromagnet in equilibrium. Within this 
spin model we investigated the dynamical reduction of 
quantum physics to classical behavior via spontaneous 
symmetry breaking. The goal of this paper is to present 
a self-contained, detailed description and explanation of 
the decoherence process caused by spontaneous symme- 
try breaking. 

The fact that spontaneous symmetry breaking can lead 
to decoherence on a time scale tspon might come as a 
surprise. For a macroscopic body at room temperature, 
h/{kBT) is of order 10~^^ seconds, which is quite a short 
time. However, multiplying with Avogadro's number 
N ~ 10^"*, yields tspon — 10^° seconds, corresponding 
to a couple of centuries. Given all other sources of deco- 
herence for such a large macroscopic body, this is surely 
not a relevant timescale. However, the mesoscopic quan- 
tum qubits of contemporary interest are typically much 
smaller and the intrinsic coherence time might be reached 
in the near future. The counterintuitive feature of this 
intrinsic decoherence mechanism linked to classical equi- 
librium is that it starts to matter when systems become 
small. 

The many-particle qubits that motivate us to study 
decoherence due to spontaneous symmetry breaking are 
realized in a number of mesoscopic solid state systems. 
For instance, by engineering aluminum on a sub-micron 
length scale, superconducting flux qubits and Cooper 



pair boxes can be manufactured. The flux qubit is a 
Josephson device that can be brought into a quantum 
superposition of two electrical currents: a left and a right 
circulating current Typically this current is car- 

ried by iV 10^ Cooper pairs. A Cooper pair box on 
the other hand is a superconducting island, containing 
N ~ 10* electrons, which can be brought in superpo- 
sition of two states with a different number of Cooper 
pairs 0, 1^. Magnetic many-particle qubits are for in- 
stance realized in molecular nanomagnets. Molecules 
with large magnetic moments can be brought into a su- 
perposition of directions of magnetization. A well stud- 
ied example is Mn^^ acetate, a molecule that contains 
12 manganese atoms, coupled together to form a total 
spin of S = 10. The molecule can be brought into a 
superposition of states with = +10 and = —10 
and coherent Rabi-oscillations of the magnetization are 
observed 0. An even larger molecule is ferritin, that 
contains about 4500 Fe'^^ ions . If the total magnetic 
moment of a ferritin molecule is brought into a coher- 
ent superposition, this corresponds to a superposition of 
N ~ 10^ spins. For these mesoscopic superconducting 
and magnetic qubits the limit in coherence due to spon- 
taneous symmetry breaking is relevant. 

This paper is organized as follows. We first introduce 
the general notion of spontaneous symmetry breaking 
and the thin spectrum. We illustrate these concepts by 
the elementary example of a harmonic crystal that breaks 
translational symmetry. We then show how the presence 
of a thin spectrum can cause decoherence. Again we 
clarify the concept by describing a theoretical setup in 
which we use the crystal as a qubit. After that we switch 
to the more involved example of the (non-abelian) anti- 
ferromagnet. We first recapitulate the properties of the 
Lieb-Mattis model and determine its symmetric and 
symmetry-broken eigenstates. Subsequently we will use 
the model to determine the effect of spontaneous sym- 
metry breaking on quantum coherence. For this purpose 
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we devise a generic gedanken experiment by coupling a 
two spin singlet state to the macroscopic Lieb-Mattis an- 
tiferromagnet. We calculate the exact time evolution of 
this many-body quantum superposition. Our main re- 
sult on decoherence follows from the evaluation of the 
reduced density matrix of the superposition by tracing 
out the thin states. Finally we will then discuss what 
it reveals regarding spontaneous symmetry breaking in 
more general Heisenberg-like spin Hamiltonians. 



Spontaneous Symmetry Breaking 

Due to the homogeneity of space, the laws of nature 
possess translational invariance. This invariance implies 
the classical law of conservation of momentum. Like- 
wise, space being isotropic enforces rotational symmetry 
on the laws of physics, implying conservation of angular 
momentum. In quantum mechanics the power of symme- 
try is even greater: the translational and rotational in- 
variance of the laws of nature taken at face value, should 
imply that any physical quantum object that obeys these 
laws, has translational and rotational symmetry. How- 
ever, daily experience shows that this conclusion is non- 
sensical. If the universe around us and everything in 
it would be translationally and rotationally invariant, it 
would look the same at all places and in all directions: 
we would be surrounded by a " quantum mist" , while hu- 
man observers should be dissolved in this quantum fog 
as well. However, in the real world translational and 
rotational symmetry are manifestly broken. The funda- 
mental difference between quantum and classical physics 
lies in the role of symmetry. Dealing with an exact quan- 
tum mechanical eigenstate, all configurations equivalent 
by symmetry should have exactly the same status in 
principle, while in a classical state one of them can be 
singled out. In the example above, given that space is 
translationally invariant, a quantum object should be in 
an eigenstate of total momentum, being spread out with 
equal probability over all of space. In the classical limit 
however it takes a definite position. The explanation of 
this 'spontaneous symmetry breaking' as a ramification 
of the singular nature of the thermodynamic limit is one 
of the central achievements of quantum condensed mat- 
ter physics j^. One imagines a symmetry breaking 'order 
parameter field' B (e.g., a potential singling out a specific 
position in space). Upon sending B to zero before taking 
the thermodynamic limit {N — > oo) one finds the exact 
quantum groundstate respecting the symmetry. How- 
ever, taking the opposite order of limits one finds that 
the classical state becomes fact. Although the concept 
of spontaneous symmetry breaking was originally intro- 
duced in the context of quantum magnetism in solid state 
physics 0, spontaneous symmetry breaking is a general 
phenomenon, that is just as relevant in other fields, in- 
cluding elementary particle physics and cosmology [lOj . 



Let us first consider how spontaneous symmetry break- 
ing arises in a crystalline lattice to continue in the next 
section with antiferromagnets. Consider the textbook 
example of a harmonic crystal, with the Hamiltonian 



H = 
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where j labels all N atoms in the lattice, which have 
mass m, momentum and position Xj. The harmonic 
potential between neighboring atoms is parameterized by 
k; it turns out that the results on spontaneous symmetry 
breaking that follow are equally valid for an-harmonic 
potentials. Let us first identify the collective dynamics 
which describe the spontaneous symmetry breaking of 
this short-ranged microscopic Hamiltonian. 

In the standard treatment of the quantum crystal one 
begins by introducing new coordinates, which are the 
displacements of atoms from their equilibrium position. 
Then, after a Fourier transform the eigenstates of this 
Hamiltonian are easily found. We take a slightly longer 
route by introducing bosonic phonon operators from the 
very beginning and diagonalizing the quadratic part of 
the Hamiltonian by performing a Bogoliubov transforma- 
tion at the end. In doing so we do not have to introduce 
any equilibrium position of the atoms. Instead we can 
keep track of the center of mass motion of the crystal as 
a whole, and this brings to the fore the thin spectrum in 
a natural manner. Moreover we can use the exact same 
procedure in the next section to find the collective order 
parameter dynamics for antiferromagnets. 

The momentum and position operators are expressed 
in terms of bosonic operators as 



so that the commutation relation [a;j,Pj'] = ifiSjji is ful- 
filled. We choose = \/2mK so that the Hamiltonian 
reduces to 



and after a Fourier transformation 



H ^ h 



K ^ — ^ 



Akblb 



where Ak — 2 — cos(fca), Bk = — cos (fca) and a is the 
lattice constant. This Hamiltonian is still not diagonal, 
since the terms ^^^L^, and create and annihilate 

two bosons at the same time. We get rid of these terms by 
a Bogoliubov transformation (see appendix). After this 
the Hamiltonian in terms of transformed boson operators 
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= cosh(wfe)6_j, + sinh{uk)bl is 



H 



2sin|fca/2| (plPk+n 



kl-'k ■ 2 



COS (ka) 



nk 



(4) 



since cos k — J^^ dk cos fc = 0. 

This result seems to coincide with the textbook Hamil- 
tonian which we would have obtained if we had fol- 
lowed the conventional route of Fourier transforming the 
Hamiltonian for the displacements, and then quantizing 
it. However, the Bogoliubov transformation has the ad- 
vantage that it brings to the fore a rather subtle point. 
When fc the excitation energy Wfe — > and the two 
parameters in the Bogoliubov transformation diverge: 
sinh(Mi.) — > CXI and cosh(ufc) —^ oo. Precisely at fc = 
the canonical transformation is no longer well defined. 
We therefore should investigate the bosonic terms in the 
Hamiltonian with k ^ separately. This zero momen- 
tum part of the Hamiltonian describes the obvious fact 
that the crystal as a whole carries a kinetic energy asso- 
ciated with the combined mass of all of its constituents, 
and is given by 
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1 
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i?k=o = + constant. 



(5) 



(6) 



where ptot = Pi = Pk=o is the total momen- 
tum of the entire system, or equivalently, its center of 
mass momentum. When iV is large, this Hamiltonian 
has states that are very low in energy. These states in 
fact form the thin spectrum of the harmonic crystal. We 
call this part of the spectrum thin because it contains so 
few states of such low energy that its contribution to the 
free energy in the thermodynamic limit completely dis- 
appears (see appendix). In turn, this implies that these 
thin spectrum states do not contribute to any thermody- 
namically measurable quantities such as for instance the 
specific heat of the crystal. Their effect on the properties 
of the crystal is thus increasingly subtle, but its existence 
can nonetheless have profound consequences. In classi- 
cal systems this thin spectrum is absent: it is quantum 
mechanics that generates it. About a decade ago the 
deep meaning of the thin spectrum for interacting quan- 
tum systems became clear and conseque ntly its explicit 
mathematical structure was determined Il2l| . 



The groundstate of the Hamiltonian at k = 0, which 
governs the collective behavior of the crystal as a whole, 
obviously has total momentum zero. It thus has no un- 
certainty in total momentum and maximum uncertainty 
in total position: translational symmetry is unbroken. 
Symmetry breaking can occur if we add to the Hamilto- 
nian of equation © a symmetry breaking field of the 
form _Bxjo(/2, where the center of mass coordinate is 
xtot = X^j^i- This yields a harmonic oscillator equa- 
tion for the collective position coordinate. Its well known 
groundstate wavefunction is 



i)o[xto 



This state describes a wave-packet for the center of mass 
coordinate in real space, which of course corresponds to 
an equivalent superposition of total momentum states: 
the symmetry breaking field B couples the different thin 
spectrum states of the crystal. For a vanishing symme- 
try breaking field B and finite number of atoms we have 
LoN and the collective coordinate is completely de- 
localized, as before: ^'0(2^*0*) = const. But taking the 
thermodynamic limit {N — > 00) m presence of a finite 
symmetry breaking field gives uN —> 00 and the cen- 
ter of mass position becomes completely localized in the 
center of the potential well {ipoixtot) = <5xtot,o), even if 
at the end the symmetry breaking field is sent to zero. 
As we already pointed out, such a singular limit charac- 
terizes spontaneous symmetry breaking; in this particu- 
lar case the translational symmetry of the crystal as a 
whole is spontaneously broken. The occurrence of a thin 
spectrum which consists of the states associated with the 
quantum mechanics of the macroscopic body as a whole 
is a universal notion. Whenever a system exhibits a con- 
tinuous symmetry which is broken in the classically re- 
alized, macroscopic state, then consequently there must 
be a spectrum of states associated with the symmetry- 
restoring fiuctuations of the orderparameter as a whole. 
The universality of this notion is guaranteed by the Gold- 
stone protection which holds for the low energy excita- 
tions at /c — !■ in any system with a broken continuous 
symmetry. The smallness of the energy spacing within 
the thin spectrum warrants the orderparameter dynam- 
ics of macroscopic bodies to take place on a time scale 
much larger than anything observable. 



Decoherence 

To study the effect of the thin spectrum on the coher- 
ence of many particle qubits, let us first investigate the 
dynamics of such a qubit in the most general terms. Con- 
sider a many particle system that is large enough to dis- 
play a spontaneously broken continuous symmetry, but 
small enough to be used as a qubit. This qubit will then 
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have a thin spectrum which we can label by the quantum 
number n. At the same time the system must have two 
accessible quantum states that can be used as the qubit 
states, and which can be labeled by the quantum num- 
ber m. Because we have no experimental control over 
the thin spectrum states, we will have to start out the 
experiment with a thermal mixture of those states: 



Pt<ta 



l^e-'^^o"|o,„) (0,n| 



(8) 



where pt is the density matrix, is the energy of the 
state \m,n) and where, by definition, the partition func- 
tion is Z = E„e"'^'^" and (3'^ = ksT. To begin us- 
ing this qubit in a quantum computation we will typi- 
cally have to prepare it in some coherent superposition 
of the states in the two level system. To do this we ap- 
ply a rotation that takes the state |0, n) into the state 
•y/l/2 (|0, n) -I- for all values of n. The resulting 

density matrix then is given by 



Pt = to 



^^e-^^?(|0,n)(0,n| + |0,n)(l,n| 



2Z 

n 

+ |l,n)(0,n| + |l,n)(l,n|). 



(9) 



If we know the Hamiltonian H which governs the dynam- 
ics of the qubit, then we can follow the time evolution of 
this density operator by applying the time evolution op- 
erator U\m,n) = e~'fe^'^*~*''^|m, n) = e~'^-^™(*~*o)|m, n). 
We then find for the density matrix at t > to 

Pt>to = Upt^toU^ 
n 

+ [e-^(^?-^")(*-*°)|0,n)(l,n| +/I.C.]) ,(10) 



where the off-diagonal matrix element is defined as 



Pt>to — 



(13) 



If this off-diagonal matrix element vanishes at some 
time, then the qubit will have decohered at that time, 
due to the presence of the thin spectrum. In general 
^Ethin = Eq — i?" will not be zero, and this shift in 
energy corresponds to a phase shift of the thin spectrum 
states. These phases will typically interfere destructively, 
lowering and leading to dephasing and decoherence. 
The timescale for this decoherence process is set by the 
inverse of the involved energy scale, and will therefore be 
proportional to h/AEthin- 

For this dephasing to occur however, it is necessary 
that a finite number of thin spectrum states participates 
in the dynamics of decoherence. How many states do 
contribute to the process is governed by the Boltzmann 
factor e~^^o ^ which exponentially suppresses states of 
energy higher than ~ ksT/EtMn with EtMn the typi- 
cal level spacing of the thin spectrum states. Putting 
these arguments together, one finds that the character- 
istic timescale on which pf-^to ^^^^ vanish should be pro- 
portional to 
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(14) 



In the following sections we will calculate tspon ex- 
plicitly for a number of realizations of the many par- 
ticle qubit. We will see that in the generic situation 
AEthin oc Ethin/N so that we find 



t 



Nh 



(15) 



which is our main result. 



where h.c. denotes the hermitian conjugate of its preced- 
ing term. 

Experimentally the thin spectrum is as good as unob- 
servable because of its extremely low energy and its van- 
ishing thermodynamic weight (see appendix). We there- 
fore have to trace these states out of the density matrix. 
This will yield a reduced (observable) density matrix, de- 
fined by 



(11) 



where the trace is over thin spectrum states labeled by 
i and (j|to, n) = \m)5j^n- Performing the trace, we find 
the following reduced density matrix in the basis of the 
states |m) 



red 
Pt>to 



OD 
Pt>to 



(12) 



THE CRYSTAL AS A QUBIT 

As a first example of the influence of the thin spectrum 
on coherence, let us try to employ the harmonic crystal 
discussed in the introduction as a qubit. In order to do 
so we will have to define a set of two states that are 
to be used as the calculational states of the qubit. A 
simple choice for such a set could be to use the presence 
or absence of an interstitial excitation. This leads to 
the definition of the state |m = 0) describing the crystal 
with N atoms, and the state |m = 1) which has one extra 
interstitial atom, and is described by the same model, but 
with A^ -|- 1 atoms in the lattice. The thin spectrum is 
exactly as described in ((HJ, so that the energy can be 
defined as 



e: 



2M {N + m) 



(16) 
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where M is the mass of an atom, n labels states with 
different total momentum (which make up the thin spec- 
trum), and /i is the chemical potential associated with 
adding an extra atom to the lattice. 

We are now in the position to simply substitute this in- 
formation into the general expression for the off-diagonal 
matrix element of the reduced density matrix H13() , yield- 



ing 



Pt>to 



h 2M \ N 



(w- jvTt)(*-*o) 



(17) 



The constant phase factor e~"s:'^(*~*o) does not contribute 
to the decoherence process, but the terms depending 
on n introduce phase shifts into the dynamics of the 
system, which lead to the disappearance of over 
time. Upon introduction of Ethin = 1/{2MN) and 
AEtMn = 1/{2MN{N + 1)) ~ 1/{2MN^), a straight- 
forward evaluation of the sum over thin spectrum states 



yields tspon , defined as the half time for 



\Pt>to I 



2TTh Et 



hin 



^spon 



UbT AEthin 

2Tih 



N 



(18) 



By using the crystal as a qubit in this way we have 
assumed that we can just ignore the symmetry breaking 
field as soon as the crystal has been localized in space at 
some time in the past. In general this may not be true, 
because the thermodynamic limit and the limit of dis- 
appearing localization field do not commute. We should 
therefore also consider the situation in the presence of a 
small but finite symmetry breaking field B. In that case 
the energies of the system will be given by 



El 



B 



2M [N + m) 



/im. 



(19) 



which will again lead to a phase factor which is constant 
in n, and a sum over phases which can be written as 
multiples of Ethm = \/B/2M {N + m) and AEtMn ~ 
Ethin/N. The summation over thin spectrum states will 
thus again yield the coherence time tspon — ^j~t- 



Goldstone Modes 

The interstitial excitation used in the previous section 
to characterize the qubit state of the harmonic crystal is 
a very rough excitation to use for that purpose. The 
extra atom in the crystal will immediately affect the 
whole lattice structure, and thus couple to many phonon- 
excitations. On top of that it also automatically increases 
the mass of the crystal as a whole, and thereby influ- 
ences the thin spectrum directly. Because of this direct 



alteration of the crystal properties, serious decoherence 
effects are to be expected. In constructing a qubit using 
the quantum crystal it is therefore better to look for a 
more 'silent' excitation. These silent excitations are nat- 
urally found in the long wavelength Goldstone modes of 
the crystal, i.e. the low energy phonons. 

To see what is the effect on the thin spectrum of having 
phonons present, let's consider the general anharmonic 
term 



(20) 



which supplements the Hamiltonian and introduces 
phonon-phonon interactions in its dynamics. The Fourier 
transform of Han is 



(21) 



with a the lattice constant. From this expression it is im- 
mediately clear that for all values of z the k = part of 
Han will vanish. This is of course just a manifestation of 
the Goldstone protection of the long wavelength phonon 
modes, but it implies that the presence of phonons can 
not alter the form of the thin spectrum part of the Hamil- 
tonian given by © . The addition of a symmetry breaking 
field will not modify that conclusion, because the ener- 
gies of the collective model in presence of a localization 
field only depend on the total mass of the crystal and the 
size of the symmetry breaking field. The phonons do not 
influence the mass of the quantum crystal, and thus we 
find that AEthin = both for zero and finite symmetry 
breaking field B. 

Because the phonons do not influence the thin spec- 
trum at all, there will also be no induced decoherence 
of a superposition of different phonon states by the thin 
spectrum. In this case thus t^pon — > cxd. Of course the 
superposition of phonon states will be influenced by all 
the other sources of decoherence that are present in the 
crystal, but it is protected from the decohering effects 
of the thin spectrum by the fact that phonons do not 
influence the properties of the crystal as a whole. 

The fact that the Goldstone modes of the crystal sup- 
ply an excitation that is protected from all decohering ef- 
fects of the thin spectrum is due to the fact that we only 
considered translation symmetry breaking of the crystal 
as a whole. Within it, there is still translation invariance 
of the atoms. This pathology that is due to the periodic 
boundary conditations that we implicitely assumed. A 
coupling between the thin spectrum and the Goldstone 
modes can appear as soon as this internal symmetry is 
broken as well -which is the generic situation. In fact 
we will show in the next section that even in the most 
silent non-Abelian model that one can think of, i.e. the 
antiferromagnetic Lieb-Mattis model, this pathology will 
disappear. 
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LIEB-MATTIS MODEL 

Let us now turn to the discussion of the antiferromag- 
netic Lieb-Mattis model. The reason for considering the 
rather particular, long ranged Lieb-Mattis model is that 
for a broad class of Heisenberg models with short-ranged 
interactions it constitutes the effective Hamiltonian for 
the thin states. Similar collective models underly the 
breaking of other continuous symmetries, such as for in- 
stance gauge symmetry in a superconductor. In that case 
the collective Hamiltonian turns out to be very similar 
to the Lieb-Mattis Hamiltonian as far as the structure of 
the thin spectrum and the composition of the wavefunc- 
tion of the symmetry broken state are concerned. To ex- 
plicitly show how the Lieb-Mattis model is arises from a 
Heisenberg model, let us consider an antiferromagnet on 
a bipartite lattice with isotropic nearest neighbor interac- 
tions between quantum spins of size a. Its Hamiltonian 
is 

i7 = J^SA+5, (22) 

where i labels all the spins on the A sublattice, and the 
(5 are the vectors connecting site i to its neighbors on 
sublattice B. The generalization to other types of in- 
teractions and even othe r ty pes of lattices is straight- 
forward HI H E E El El- The magnon spec- 
trum of this Hamiltonian can be found within linear spin 
wave theory. One approximates the spin operators with 
Holstein-Primakoff bosons as 

SleA ^ cr - a\a.^ , Sl^Q h\h^ ~ cr, 

Sr^^ V2^a\, Sr^^V2^b,. (23) 

To quadratic order in the boson operators the Hamilto- 
nian becomes, after a Fourier transformation 

H^S"^ =^-JNza' + JzaJ2[iala^ + blb^) 

k 

+ 7k(4&^-k + «k^-k)). (24) 

where z is the coordination number of the lattice, N the 
number of lattice sites and 7^, = ^^ge^^^. The last 
two terms in this expression can be diagonalized by a 
Bogoliubov transformation (see appendix). Again the 
important point is that the Bogoliubov transformation is 
singular at fc = and A: = tt, as in both cases 7^ — > 1. We 
therefore treat these two fc-points separately. Turning 
back to the notation in terms of spins, using that the 
Fourier transform of our Hamiltonian is 

H = JY.^A-^-^ (25) 

k 



and 

Sk=7r = ^=(Sa — Sb) 

V N 

we find that the singular parts of the spectrum reduce 
exactly to the Lieb-Mattis Hamiltonian 

H = + J ^kSk • S_k 

k7tO,ir 

9 T T 
HZ7 = f • Sb = ^(S^ - - S|), (26) 

where and are the total spins of each sublattice, 
and S is the total spin of the system. From here on we 
will focus entirely on this collective Hamiltonian, as it is 
the only part of the Heisenberg-like Hamiltonians that is 
relevant for the spontaneous symmetry breaking of the 
antiferromagnet as a whole. Notice that the internal or- 
dering of the individual spins within the antiferromagnet 
can be destroyed by fluctuations of finite wavelength that 
we do not consider in this collective, long ranged model. 
We assign to the Hamiltonian H'^j^ the superscript sym 
because this Hamiltonian, as we will show below, de- 
scribes the symmetric (symmetry unbroken) state of the 
antiferromagnet . 

In H^^^ each spin on the A sublattice interacts with 
all spins on the B sublattice and vice versa, thus creating 
infinite range interactions. The energies of the Hamilto- 
nian are trivially identified as [S{S +1)- Sa{Sa + 1) - 
SsiSs + 1)] and the corresponding eigenfunctions are la- 
beled by their quantum numbers \Sa, Sb, M). Here 
the z-component of the total spin S is denoted by M. 
Clearly the ground state of Hf^^ is a singlet of total 
spin: the state with lowest energy has S — 0. In fact 
there is an exact proof that the groundstate of any finite 
spin system of this sort is a total spin zero (5* = 0) sin- 
glet [l9|. Notice that all states which differ only in M 
are degenerate. For simplicity (and without loss of gen- 
erality) we henceforth take the quantum number M to 
be zero 0. 

The groundstate singlet \Na/2, Na/2, 0, 0), with both 
Sa and Sb maximal and S = M = Q is separated by 
energies of order J /N from states with higher S. The 
set of these extremely low energy states that only dif- 
fer in their total spin quantum number forms the thin 
spectrum 0, 0, . Since is contained in as its 
k = and k = tt components, and since the thin spec- 
trum of the Lieb-Mattis model is formed by the k = 
component, exactly the same thin spectrum must govern 
the collective dynamics of other antiferromagnets with 
short-range interactions 0, 0, 0| . 

There are also excitations in (|26|l that can be created 
by lowering Sa or Sb- This costs an energy of order J, 
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and it can easily be shown that these excitations corre- 
spond to the elementary excitations, the magnons, of the 
Lieb-Mattis system |20j. Because of the extremely long 
ranged interactions the magnons are gapped and disper- 
sionless. 



The Lieb-Mattis Hamiltonian with symmetry 
breaking field 

Having defined the Lieb-Mattis model in its symmet- 
ric form, we now review how to explicitly break its non- 
Abclian SU(2) spin rotation symmetry. We will show 
that in the thermodynamic limit the symmetry breaking 
occurs spontaneously 0| . Since the groundstatc of H'^j^/ 
is a singlet of total spin, this state is orthogonal to the 
Neel state, which is the ground state of a classical antifer- 
romagnet. We should stress here that there is a marked 
difference between ferro- and antiferromagnets. Even if 
spontaneous symmetry breaking is very often discussed 
with the example of a ferromagnet at hand, the spon- 
taneous symmetry breaking in a ferromagnet is not the 



generic situation for an interacting quantum system. The 
reason is that the total magnetization (pointing along, 
e.g., the z-axis), which is the order parameter of a ferro- 
magnet, commutes with the Hamiltonian: it is nothing 
but the projection of the total spin along that axis, 3^^^. 
So the situation arises that the orderparameter is a con- 
stant of motion, which is a pathology of the ferromagnet. 
This same pathology leads to the absence of an interest- 
ing thin spectrum, because in the ferromagnet states with 
different S^^^ are strictly degenerate. Quantum systems 
in general, however, have non-trivial thin spectra. 

Refocusing on antiferromagnets, we need to proof that 
the Neel state is a stable groundstate in the thermody- 
namic limit. In order to do so, an explicit symmetry 
breaking field B is introduced : 



H 



LM 



TTSym 



B{S\-S%). 



(27) 



Clearly the symmetry breaking field induces a finite sub- 
lattice magnetization. The field couples the different to- 
tal spin states of the thin spectrum by the matrix ele- 
ments 



J 



{Sa, Sb, S, M\S\ - S%\S'a, S'b,S', M') = 5sa,s' Jss,s' s5m,m' [fs+i5s,s'-i + 9s5s,s' + fs5s,S'+i] (28) 



where 



, _ \S''-{SA-SBf][{SA + SB + lf~S^][S^-M^] , _{Sa-Sb){Sa + Sb + 1)M 



{2S +1){2S -1)S^ 



S{S+1) 



These matrix elements are found by performing a rather 
tedious sum over Clebsch Gordon coefficients in the fol- 
lowing expression (see also the appendix): 



{Sa, Sb, S, M\S% — Sb\S'a, S'b,S', M') — 

X! \^s'a!sbMaM-Ma'~'Sa^SbMa,M-Ma^'^^'^A ~ 
Ma 



SSa^S' a^Sb,S' B^M,M' 



(30) 



The spectrum of cigenstatcs \n) in the presence of a 
symmetry breaking field can now be found by expand- 
ing these states in the basis of total spin states: \n) — 
12s ^sl*^) i^'^^ clarity of notation we suppress the depen- 
dency of Ug and other variables on the quantum numbers 
Sa, Sb and M). In this basis, Schrodinger's equation be- 
comes 



HLM\n)^ES\n) ^ 



JS(5+1) „ 
N 



E 



LM "S 



where Ef^^ is the groundstate energy of H^'j^ and E^ 
is the energy of eigenstate \n) of the symmetry broken 
Hamiltonian -its thin spectrum. Here we restricted our- 
selves to the zero-magnon subspace, where Sa = Sb — 
Na/2 (hence the subscript in Eq). The generaliza- 
tion to systems with a finite number of magnons will be 
straightforward. 

In the continuum limit where N is large and <^ S <^ 
N, the matrix elements due to the symmetry breaking 
field simplify considerably. It is easy to show that in this 



Ncr/2. (32) 




We will see shortly that only the first « vTV total spin 
states contribute to the groundstate wavefunction, so 
that an expansion in S/N is justified. Notice that when 
the sublattice spin Sa is reduced by one, i.e. when there 
is a spin- wave present, the matrix element fs is reduced: 
fs ~ =fsi^-w^), for large N. This reflects 

(31) the fact that a magnon reduces the Neel order parame- 





FIG. 1: Wavefunctions of the thin-spectrum state in presence 
of a symmetry breaking field, in the continuum Umit. The 
boundary condition S > implies that of the harmonic oscil- 
lator solutions (left) only the odd ones are allowed (right), as 
these have a node at the origin. 



ter (the staggered magnetization) by unity. This effect is 
small, but turns out to be essential when we shall con- 
sider the quantum coherence of magnons: dephasing will 
occur because magnons give rise to a subtle change in the 
level splitting of the thin spectrum. This change in level 
splitting turns out to be inversely proportional to N, the 
total number of spins in the antiferromagnet. 

In the continuum limit the Schrodinger's equation 
reduces to 



1 d 



1 



(33) 



where again we have used S N. In this equa- 
tion we introduced i 



^Vli and ly,, = 



BNcr 



1. 



Obviously this is the differential equation of a harmonic 
oscillator. The eigenstates thus are well known and 
the corresponding eigenvalues are Vn = {n + 1/2)^^1 so 
that 



1 1 F 

2 I ^thin^ 



(34) 



where the quantum of energy for the states labeled by n 
is Ethin = \/2aJB. For the harmonic oscillator n is a 
non-negative integer. However, in the present situation 
we have to meet the boundary condition that 5* > or, 
equivalently, that = if S* < 0. So has to vanish 
at the origin. This boundary condition is trivially met by 
eigenf unctions that are odd and have a node at S* = 0, see 
Fig.n Thus solutions to the Schrodinger's equation H33|) 
are harmonic oscillator eigenfunctions of order n, where 
n is an odd positive integer. In Fig. |21 the groundstate 
wavefunction in the continuum limit is compared with 
the exact wavefunction for large N . It makes clear that 
the continuum approximation is very good one. 

Let us consider the energy spectrum in Eq. 1)34(1 in 
more detail. Clearly if B is zero we recover the ground- 
state energy Ef^^ of the symmetric case that we dis- 
cussed before. However, if there is a finite staggered 
field B, there is a gain in groundstate energy propor- 
tional to BN , which reveals that the energy spectrum in 
Eq. H34|) is the one of a Neel state. The same conclusion 
is reached by directly calculating the orderparameter (see 




->S/N 



FIG. 2: Comparison of the exact symmetry broken Neel wave- 
function (for N = 500 spins) and the Neel wavefunction in 
the continuum limit (a harmonic oscillator eigenstate). The 
overlap of the Neel state with the different total spin states 
is shown as a function of the total spin quantum number. 
The parameters are J = 1 and the symmetry breaking field 
is B = 1/10. The wavefunctions are rescaled such that the 
maximum of the harminic oscillator wavefunction is unity. 



appendix). The result is shown in Fig.|2| Apparently, for 
the symmetry broken Neel state to be stable, the symme- 
try breaking field can be exceedingly small, as long as N 
is large enough. In other words: in the thermodynamic 
limit the spin rotation symmetry of H^^^ can be sponta- 
neously broken by an infinitesimal field B. Putting it in 
a more formal manner: spontaneous symmetry breaking 
gives rise to the singular limit 



lim lim 



lim lim 



^B 

Na 

'^B 

Na 



and 
1. 



(35) 



This in fact defines spontaneous symmetry breaking, just 
as it did in the case of the quantum crystal. That for the 
Lieb-Mattis Hamiltonian this limit is singular is directly 
clear from Fig. |3| 

In the symmetry broken Neel state the excitations la- 
beled by n act as a new thin spectrum with excitation 
energies that are multiples of Ethin — V2aJB. The 
magnon excitation energy is still of order J. 

We now repeat the analysis above for a Neel state with 
m magnons (by setting Sa = Sb = Na/2 — m/2 and 
using /™ instead of /s). In this case the energy spectrum 
becomes 



E'^^E^ + m{2aJ + B) 



2Na 



{n+l/2)EtMn. (36) 



see Fig. 01 Note that as we stated before, there is a subtle 
effect of the magnons on the thin states: m magnons 
cause a change in energy of the thin spectrum of the 
order of m/N . This effect will turn out to be essential 
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Orderparameter versus symmetry breaking field 



Semi-Classical Qubit 



Z 

. 

A 



00 
V 



t<t„ 



exact 

continuum 



-6/N" 2/N' 6/N^ 

-> B/J 



FIG. 3: Orderparameter as a function of symmetry breaking 
field. The exact result for A'^ = 100 spins, the continuum 
expression for the orderparameter is derived in the appendix. 



Energy level diagram 
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FIG. 4: Energy level scheme with the zero and two magnon 
states, each with its tower of thin spectrum states. The level 
spacing in the thin spectrum is Ethin, magnons live on an 
energy scale J. 



for the decoherence mechanism that we will discuss in 
the following section. 

Physically the change in energy of the thin states due 
to the presence of a magnon can easily be understood. 
If there are m magnons present in the antiferromagnet, 
then the order parameter of the total system is reduced 
by TO. Since the thin spectrum describes the global exci- 
tations of the order parameter, its energy is proportional 
to the order parameter itself. The ratio of the Neel or- 
der parameter of the excited state with to magnons and 
the one of the groundstate with a fully developed order 
parameter is {N — m)IN . Therefore, when there are ni 
magnons present, the relative change of the orderparam- 
eter is to/TV and the change in energy of the thin states is 
therefore of the order niEthin/N which explains the last 
term in expression (|35|l . 
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\0)-\2) 


Rabi Oscillations 
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FIG. 5: Semi-classical time evolution of a two spin qubit that 
at t = to starts interacting with a A^-spin Lieb-Mattis mag- 
net, thus forming ai t > to a. many-body qubit made out of 
N + 2 spin. Quantum coherence is preserved at all times, 
since in this semi-classical approximation the thin spectrum 
is neglected. 



Preparing the many-spin qubit 

Using the many-body Lieb-Mattis model with TV spins 
and (7 — 1/2, we now study the coherence of the anti- 
ferromagnet when it is used as a qubit. Again there are 
many ways in which one can define a two-level system 
to be used as the qubit states. The best possible choice 
in this case is provided by the gapped and dispersionless 
magnons: we use as a qubit (or cat-state) the super- 
position of a perfectly ordered antiferromagnet and the 
state of the antiferromagnet with one magnon on each 
sublattice. Due to the long-range nature of the interac- 
tion in the Lieb-Mattis model the gapped magnons them- 
selves are not damped and as such do not decay or deco- 
here. Also, in analogy to the quantum crystal we expect 
the magnons or Goldstone modes of the antiferromagnet 
to influence the thin spectrum as little as possible. A 
magnon has an energy J, which we assume is the energy 
scale that is available to the (thought-)experimentalist to 
prepare, manipulate and read out the qubit |2ll |. 

With the exact expressions for all eigenstates and en- 
ergies of both the symmetric and the symmetry broken 
Hamiltonian at hand, we are in the position to set up 
the initial state for our many-particle qubit. Instead of 
simply assuming that we are in a previously prepared 
superposition of states with zero and two magnons, we 
will explicitly construct this initial state. This can be 
done by coupling at time t — to a two spin singlet to the 
symmetry broken A^-spin Lieb-Mattis system, see Fig. 
From this the desired superposition results. So for times 
t < to, the Lieb-Mattis antiferromagnet is completely 
decoupled from the two spin singlet and the total wave- 
function is thus the direct product of the wavefunctions 
of the iV-spin magnet and the two-spin singlet state: 



\A<to) = |0,n) [singlet) 



(37) 



Here we denote the Neel state with to magnons and n 
thin spectrum excitations by |TO,n). The state \singlet) 



is 



V2 



[I iiT2) ^ I Tii2)]- Upon instantaneous inclusion at 



t = to oi the two spin state in the Lieb-Mattis lattice. 
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the groundstate of the decoupled system at i < io can be 
expressed in terms of the eigenstates of the N + 2 spin 

I 



system at t = Iq. The exact groundstate wavefunction is 
then given by the following formidable expression 



N-l 



N-1 



IV't<to) = J2^s \S, 0) (E) \singlet) = 



s=o 



^N/4:,N/4,Ma,-Ma 1/2,1/2, Ml, -Ml N/4,1/2, Ma, Ml 

s=o 

^Sb2, -Ma-Mi 



^0,0 



^Sai,Ma+Mi 







^N/4,1/2.-Ma]-Mi^Sai,Sb2,Ma + Mi,-Ma-Mi^S,St\Sai, Sb2, S, 0), (38) 



where we sum the Clebsch Gordon coefficients over Ma, 
Ml and over the total spins Sai, Sb2 and S, A/B denote 
the spins on sub-lattices A and B and the spins on sites 
1 and 2 make up the singlet. With Al(B2) we denote the 
set of spins on sublattice A(B) combined with spin 1(2). 
The sums can be evaluated and we obtain 



N-l 

|V't<to) = XI 

s=o 



(W-2 + 2S)(JV + 4 + 2S) I W + 2 JV + 2 
2(JV + 2)^ I 4 ' 4 



S,0) 



I ; 2S(S + 1) I iV + 2 JV-2 
(N + 2)^ I 4 > 4 



s,o) 



2S(S + 1) I W-2 JV + 2 



S,0 



(W-2S)(JV + 2 + 2S) I JV-2 JV-2 



S,0 



(39) 



Again the equations simplify drastically in the contin- 
uum limit of large N , where as before <C 5 ^ A^. In 
this case the wavefunction of the system at t < to, ex- 
pressed in the eigenstates of the A^ -I- 2 spin system at 
t = to is 



|Vt<to) = [|0,7^) + |2,n)]/V2. 



(40) 



Here all states on the right hand side, i.e. all the thin 
spectrum states labeled by their quantum number n with 
either zero or two magnons, refer to configurations of 
N + 2 spins. 

To account for a finite temperature of our many parti- 
cle qubit, we combine initial states with different n into 
a thermal mixture before we let it interact with the two 
spin singlet. We should stress that we only consider 
temperatures that are much below the magnon energy: 
ksT ^ J so that there is no thermal occupation of the 
magnon states This implies that the order param- 
eter is not affected by the thermal fluctuations. So, all 
that we introduce is an incoherent mixture of the low ly- 
ing thin spectrum states, which all support a finite sub- 
lattice magnetization. Still, the implicit assumption is 
that the thin states are in thermal equilibrium -and it 
is an important assumption as our final result relies on 
it. In principle it can of course not be excluded that oc- 
cupation distribution of the thin spectrum states is far 
from thermal equilibrium. But as we have not a priori 
prepared the thin states in some particular way, we as- 
sume them to be thermally occupied. Then the density 



matrix at times t < is then 

Pt<to = e~^^° |0, n) (g) \singlet) {singlet] (g) (0, n\ 

n 

71 

+ \2,n)(0,n\ + \2,n){2,n\) (41) 

where, by definition, the partition function is Z = 
X;„e-'3^o" and f]-^ = fcsT. 

TIME EVOLUTION AND DECOHERENCE. 

By coupling the symmetry broken Lieb-Mattis model 
to the two spin singlet, we have created the initial state of 
our A^ + 2 spin qubit. This initial state is precisely equiv- 
alent to the initial state © of the general description, 
and we can thus follow equations IjlOl) . (|ll|l and H12|l di- 
rectly. That way we compute the exact time evolution of 
the initial state density matrix, trace away the thin spec- 
trum states which have vanishing thermodynamic weight, 
and finally define the off-diagonal element of the reduced 
density matrix as 



Pt>to — 



— V e-'3^o e-*(^o"--E2")(*-*o)/s 

7 ^ 



(42) 



We can then substitute the exact expressions for i?'" in 
this matrix element, and perform the summation. We 
find 



Pt>to ^ 



(43) 



1 — e^^^ 1 — Q—x—iT 

where x = and r = ^AEthinit - to), with Ethin = 

V JB and AEthin — Ethin/ N . We again define the co- 
herence time tspon as the half-time of | | . From the 
equation above one then finds 



27rA^;i 



(44) 



our main result. 
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Notice that just as in the case of using a quantum 
crystal with an interstitial excitation, the coherence time 
tspon in the end does not depend on any details of the 
underlying model. The fact that AEthin is proportional 
to Ethin itself removes all dependence of tspon on the 
model parameters J and B. 

Physical interpretation 

It is remarkable that the coherence time is such a uni- 
versal time-scale, independent the detailed form of the 
thin spectrum -which, after all, is determined by the pa- 
rameters J and B in the Lieb-Mattis Hamiltonian. Math- 
ematically this is due to the fact that both x and r are 
proportional to Ethin- Physically one can think of this 
universal time-scale as arising from two separate ingre- 
dients. First, the energy of a thin spectrum state |n) 
changes when magnons appear, as we pointed out above. 
The change is of the order of nEthin/N, where Ethin is 
the characteristic level spacing of the thin spectrum that 
we happen to be considering. The fact that each thin 
state shifts its energy somewhat at t > to leads to a phase 
shift of each thin state and in general these phases inter- 
fere destructively, leading to dephasing and decoherence. 
The larger nEthtn/N , the faster this dynamics. 

But from the argument above it is clear that in or- 
der for this dephasing to occur, it is necessary for a fi- 
nite number of thin states to participate in the dynamics 
of decoherence. Since temperature is finite (but always 
small compared to the magnon energy) a finite part of 
the thin spectrum is available for the dynamics. Thin 
spectrum states with an excitation energy higher than 
fc^r are suppressed exponentially due to their Boltz- 
mann weights. Therefore the maximum number of thin 
states that do contribute is roughly determined by the 
condition that ri""*^ ~ ksT / Ethin- Putting the ingredi- 
ents together, we find that the highest energy scale that 
is available to the system to decohere is approximately 
vT"^ Ethin/ N = J^^nv^- AH together, the thin spec- 
trum drops out of the equations. The time scale at which 
the dynamics take place is determined by the inverse of 
this energy scale, converted into time. From this ar- 
gument we immediately find again the coherence time 

l-spon ^ kgT ■ 

This physical picture also suggests an alternative way 
of introducing decoherence into the many particle qubit. 
Instead of raising the temperature and making an inco- 
herent superposition of more and more thin spectrum 
states, we could start out sX t < to with the Lieb- 
Mattis antiferromagnet in its (zero temperature) sym- 
metric ground state, and then instantaneously turn on 
the symmetry breaking field B sA, t = t^. At t > to 
the eigenstates are the Neel-like thin spectrum states \n). 
No magnons are created by switching on the symmetry 
breaking field. As we can expand the the states \n) in the 



basis of total spin states as |n) = ^gUglS) we can, by 
the inverse transformation, expand the total spin singlet 
state in the basis of the Neel-like thin spectrum states 
as IS" = 0) = J2n^nW)- We can now use this singlet 
state as the initial state for our qubit. This singlet state 
is a superposition of all of the different Neel-like states, 
which are separated by energies Ethin- This procedure 
thus roughly corresponds to creating a 'maximal tem- 
perature' ksT ^ N Ethin- When time evolves all of these 
states pick up different phases, which leads to decoher- 
ence when we trace over them. The coherence time due 
to this switching on of the symmetry breaking field is 
therefore tsB = ^¥=- 



Symmetric case and short-ranged models 

This raises the question what would have happened 
if we had not broken the symmetry in the Lieb-Mattis 
magnet (by introducing a finite symmetry breaking field 
B) at all. In the symmetric case S, Sa and Sb are good 
quantum numbers at all times. It is easy to see that in 
this situation the thin spectrum, determined by the quan- 
tum number 5, is independent of the "magnon" states, 
which are determined by the quantum numbers Sa and 
Sb- Since in the symmetric Lieb-Mattis Hamiltonian the 
thin spectrum does not communicate with the magnons 
and vice versa, we will find AEthin — 0, and accordingly 
no decoherence. 

The fact that S, Sa and Sb are all good quantum num- 
bers, may be regarded as a pathology of the Lieb-Mattis 
model. In fact, the model is integrable just because there 
are so many conserved quantities. In a more general, 
short ranged Heisenberg model the magnons will acquire 
a finite lifetime and it is expected that they will in gen- 
eral influence the structure of the thin spectrum, even 
if the symmetry breaking field is absent. In this sense, 
the Lieb-Mattis model can really be seen as the best case 
scenario for avoiding decoherence in non-Abelian, SU(2) 
symmetric models. Its infinitely long ranged interactions 
introduce a large energy gap for all magnons, which thus 
become extremely 'silent' excitations. On top of that the 
coupling to the collective dynamics is so subtle that it 
can only be seen because of the existence of a singular 
limit: Only because we need to always consider an in- 
finitesimal symmetry breaking field when looking at the 
thermodynamic limit do wc find decoherence at all. 



Recurrence 

Finally we notice that the off-diagonal elements of the 
density matrix, Eq. (|43|l . are periodic in time, and the 
initial density matrix recurs when Nt = 2tt or, equiva- 
lently, tree = t^^/ Ethin- Such a periodicity is required by 
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the fact that the time evolution is unitary. As the recur- 
rence time is inversely proportional to the level spacing 
of the thin spectrum, it depends on the microscopic pa- 
rameters of the model. It becomes infinitely long if the 
symmetry breaking field vanishes. In the physical limit 
the recurrence time is always much longer than the de- 
coherence time as tree /t spon = kBT/Ethin » 1- 

CONCLUSIONS 

A many-body qubit has an intrinsic limit to its max- 
imum coherence time. In this paper we have presented 
in detail the considerations and calculations that lead to 
this conclusion. 

The limit to coherence is caused by the thin spectrum. 
In quantum systems a continuous symmetry can spon- 
taneously be broken in the thermodynamic limit due to 
these states. The thin modes can be identified with the 
collective, zero momentum, excitations of the orderpa- 
rameter. In this paper we have outlined a general proce- 
dure for finding the thin spectrum states in a quantum 
systems. The states within the thin spectrum are ex- 
tremely low in energy and at the same time they are so 
few that their contribution to the thermodynamic parti- 
tion function vanishes. 

If a symmetry breaking field is introduced, the result- 
ing symmetry broken ground state is a superposition of 
(only) these thin spectrum states. The fact that the 
formation of the symmetry broken state occurs sponta- 
neously in the thermodynamic limit can then easily be 
checked by considering the non-commuting limits of dis- 
appearing field and sending the number of involved par- 
ticles to infinity. 

This has important consequences when a many-body 
quantum system is brought into a superposition of two 
different internal states. We have shown that in that case 
the thin states will in general participate in the time evo- 
lution of the full many-body system, even if their effect on 
any thermodynamic quantity vanishes. This leads to de- 
phasing and therefore decoherence when the thin states 
are integrated out. We have found that the time-scale 
corresponding to the dephasing process depends only on 
the energy scale of the thin spectrum and the energy 
shifts induced in the thin spectrum by the superposed 
initial states. Because the shifts in energy generally are 
proportional to the level spacing itself, the decoherence 
time in the end depends only on the temperature and 
size of the system, and not on the underlying details of 
the model. 

We have shown how such superpositions can be defined 
and studied in a quantum crystal and in the Lieb-Mattis 
antiferromagnet. The obvious question is to what extent 



the Lieb-Mattis qubit is representative of a general many- 
body qubit. In fact the Lieb-Mattis qubit is the best case 
scenario for the kind of many-body qubits envisaged in 
main stream quantum information theory, as its behavior 
is extremely close to semi-classical due to the presence of 
the infinite range interactions. Qubits characterized by 
short range interactions carry massless Goldstone modes 
and these will surely act as an additional heat bath lim- 
iting the coherence time. It is of course not an accident 
that the most 'silent' systems are qubits based on su- 
perconducting circuitry, which have a massive excitation 
spectrum in common with the Lieb-Mattis system. We 
have demonstrated here that even under these most fa- 
vorable circumstances quantum coherence eventually has 
to come to an end, because of the unavoidable condition 
that even the most 'silent' qubits are subtly influenced by 
their quantum origin. These effects become noticeable in 
the mesoscopic realm and we present it as a challenge to 
the experimental community to measure the maximum 
coherence time that they give rise to: tgpon ~ ^k^r ■ 



APPENDIX A: THERMODYNAMIC WEIGHT 
OF THIN SPECTRUM STATES 

It is easy to show that the contribution of the thin 
spectrum of the symmetric iV-spin Lieb-Mattis Hamilto- 
nian to the free energy density is proportional to and 
thus vanishing in the thermodynamic limit. The energy 
of a state with total spin 5* is Ethin = J/NS{S + 1) and 
its degeneracy is 25* + 1 , so that the contribution of the 
thin states to the partition function is 

N 

ZtKin = ^(25+l)e-'^^'-" 
s=o 

f-N AT 

« / (25+l)e-4^«(«+i)d5«^, (45) 

in the limit of large A''. Therefore its total contribution 
to the free energy is Fthin = —T In Zthin oc — InA^ and 
for large N its contribution to the free energy per spin 
-and the free energy density- is proportional to . 



APPENDIX B: MATRIX ELEMENTS 

By expressing the matrix elements of all components 
of the spin operators Sa and Sb in terms of Clebsch- 
Gordon coefficients, one can evaluate them by perform- 
ing the appropriate summations. The resulting matrix 
elements are: 
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(S'aS'bS'M' \S^\SaSbSM) 



Tds',S+l\l 4(4S'^-1)S'^ +OS',S 



±ds',S-l\l ^ ^ 4(4S^-1)S:^ 



2S(S+1) 
Ss',Sa^S',Sb^M',M±1 



(46) 



{S'aS'bS'M' \S\\ SaSbSM) = 

r j {S'^-{SA-SBf){{SA+SB + lf-S'^){S'^-m , {(Sa-Sb){Sa+Sb + 1)+S(S+1)) 
'>S',S+iy 4(4S'2-1)S'2 + ^S',S 2S(S+1) 



, J /(S2_(s^_Sj3)2V(s^+Sb+1)=-S2)(S2-M2) 
+ds',S-l V ^ 4(4S^-1)S^ 



Ss'y^,SA^S'g,SB^M',M 



(47) 



(S'^S'sS'M' \S^\SaSbSM) = - {S'aS'^S'M' \SI\SaSbSM) + 5s',s^{S ± M + 1){St M) 



(48) 



(S^i^J^^'M' Sa5s5M) = - (5;i5^,5'M' \S\\ SaSbSM) + M5s,,s. 



r 



(49) 



APPENDIX C: THE BOGOLIUBOV 
TRANSFORMATION 

We have used a Bogoliubov transformation to diago- 
nalize bosonic bilinear Hamiltonians of the form 



(50) 



where A^k = Ak and B^k = Bk- The relevant trans- 
formed bosons (3^f. are defined through 



h\ = cosh('u/c)/3Lj, - sinh(iifc)/3fc, 
= cosh(ufe)/3fe - sinh(ufc)/?l^. 



(51) 



The parameters Wfe obey Uk = U-k and are chosen such 
that the Hamiltonian reduces to diagonal form. This 
implies that 



cosh(2Mfe) 
sinh(2ufe) : 



A, 



and 



H = Y.^]^l-Bl U^^+^^-W (52) 
k ^ 

We can also use the exact same definition of Uk to diag- 
onalize the Hamiltonian 

H = Y,A, (4a, + blh,) + Bu [a,b_, + albl,) .(53) 



In that case the transformed bosons and Hamiltonian are 
given by 

al = cosh(Mfe)/3l;j - sinh(Mfe)a/c, 

bl = cosh(ufc)a^ - smh{uk)(3k and 

H = ^Al-Bl [alak + plpk + l)-Ak. (54) 

k 

APPENDIX D: THE ORDER PARAMETER 

Consider the symmetry broken Lieb-Mattis Hamilto- 
nian: 



H = -j^^A ■ Sb - B {Sa- Sb) 



(55) 



If the number of spins N is large, then the eigenfunctions 
|n) of this Hamiltonian are to a very good approximation 
given by the eigenfunctions of (half of) a harmonic oscil- 
lator: 



|n) = 



where \S) are the total spin cigcnstates, are the Her- 
mite polynomials, uj equals jt\/^, and n can only be 



14 



an odd integer number. Using this exact expression to 
calculate the ground state expectation value of the order 
parameter, we find: 

(^1-^I) = E44' {S'\Sl-SI,\S) 

S,S' 



= 2 XI 44-1 \ 



((7V/2 + lf 



52 52 



452-1 



(57) 



The shape of the function Ug guarantees that S <^ N, so 
that for large N the expectation value is approximately 
given by: 



N 

T 

N 

T 

N 
'2 



i: 



e 4" 



— e V iw + 0(1), 



^3/2^g-i^s^ (5-l)e-^-(^-i)'d5 
1 



1 - erf 



(58) 



which reduces to the classically expected order parameter 
in the thermodynamic limit. Note that in this expression 
it is immediately clear that the limit of vanishing sym- 
metry breaking field does not commute with the limit of 
infinitely many spins: they form a singular limit. 
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